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Abstract — In February-March 2009 NASA JPL conducted an 
airborne field campaign using the Passive Active L-band System 
(PALS) and the Ku-band Polarimetric Scatterometer (PolSCAT) 
collecting measurements of brightness temperature and near 
surface wind speeds. Flights were conducted over a region of 
expected high-speed winds in the Atlantic Ocean, for the purposes 
of algorithm development for salinity retrievals. Wind speeds 
encountered were in the range of 5 to 25 m/s during the two weeks 
deployment. The NASA-Langley GPS delay-mapping receiver 
(DMR) was also flown to collect GPS signals reflected from 
the ocean surface and generate post-correlation power vs. delay 
measurements. This data was used to estimate ocean surface 
roughness and a strong correlation with brightness temperature 
was found. Initial results suggest that reflected GPS signals, using 
small low-power instruments, will provide an additional source 
of data for correcting brightness temperature measurements for 
the purpose of sea surface salinity retrievals. 


I. Introduction 

The thermal emission of seawater in microwave frequen- 
cies is a function of the dielectric constant, which in turn 
depends upon the sea surface salinity (SSS) and sea surface 
temperature (SST). L-band brightness temperature has long 
been studied for the remote sensing of SSS[1][2] through the 
use of this principle. Brightness temperature, however, is also 
strongly dependent on the surface roughness [3] [4] and this 
dependence must be removed or corrected first, to retrieve 
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salinity measurements from microwave radiometry. In most 
cases, this correction is made using auxiliary measurements of 
surface roughness, such as those provided by an active radar 
(scatterometer). 

Bistatic radar, using signals from Global Navigation Satel- 
lite Systems (GNSS) as illumination sources, has also 
been shown to produce measurements of ocean surface 
roughness[5][6][7]. The global availability of these range- 
coded signals, transmitted on L-band frequencies (1.575 GHz) 
very near to the bands utilized for microwave radiometry, 
suggest them as potentially useful source for the roughness 
correction in microwave radiometry. This application has been 
proposed and studied theoretically[8][9][10]. In this paper, 
we will present experimental results comparing the brightness 
temperature correction using this technique (hereafter referred 
to as GNSS-R) with those from the Ku-band polarimetric 
measurements from the PolSCAT instrument. 

The GNSS-R technique has several attractive features, in 
addition to the fortuitous L-band frequency allocation for satel- 
lite navigation. Through re-utilization of existing signals with 
a passive receiver, the instrumentation would be considerably 
smaller and use less power than an active radar. Forward- 
scatter geometry results in a stronger received signal and 
the “process gain” provided through de- spreading the spread 
spectrum signal results in an increase in signal to noise ratio. 
These effects enable taking useful measurements with a low- 
gain antenna. Sensitivity to the surface roughness is present 
in the shape of the Delay-Doppler map (DDM), generated 
by cross-correlating the reflected signal with a local signal 
replica. Estimation of the probability density function for 
surface slopes can be found by inversion of the DDM, which is 
insensitive to total scattered power and thus the measurement 



does not require accurate calibration. 

Most previous work has used the dependence of brightness 
temperature on near surface wind speeds[ll][12]. In the re- 
search presented in this paper, we will use the mean square 
slope (MSS), under the assumption of an isotropic, normal 
distribution of surface slopes, as the variable of interest, and 
derive an empirical relationship between the MSS obtained 
from inverting the GNSS-R measurements and brightness 
temperature (AT*,) measured by the PALS instrument. This 
derived relationship will then be used to invert the GNSS-R 
measurements to obtain estimates of A T*,. These retrievals of 
A Tk will be compared with those obtained through applying 
the same approach to the POLSCAT wind speed measure- 
ments. Finally, we will attempt to integrate the GNSS-R MSS 
data with PolSCAT wind speed in a combined estimation of 
A T b . 

The paper is organized as follows. Section II provides a 
review of the scattering model applied to the GNSS-R signal, 
showing the relationship between the slope PDF and the shape 
of the DDM. Section III describes the experiment conducted 
February-March 2009. Section IV describes the data reduction 
and section V presents the results showing the dependence of 
MSS on A T&. Conclusions and recommendations for further 
work are given in section VI. 


II. GNSS-R Signal Model 


A sufficient model for the direct line-of- sight signal received 
from each GPS satellite takes the following form 

s(t) = A(t — r)D(t — r)p(t — r)e^ _27r ^ c+ ^ jD ^ + ^ J +n(t) (1) 


in which p(t — r) is a baseband signal with binary phase-shift 
keyed (BPSK) modulation from a pesudo-random noise (PRN) 
code and D(t) is the BPSK-modulated data message. For the 
C/A code on the LI channel (which was the signal used in all 
of the data collected in this experiment), p{) is a 1023 -chip 
PRN code with a chipping rate of 1.023 MHz, giving a 1 ms 
period in the absence of Doppler. f c is the carrier frequency 
generated at the satellite transmitter (1.575 GHz for the LI 
channel) and fn is the Doppler shift resulting from relative 
motion between the transmitter and receiver. is the phase 
which is proportional to the delay modulo 2i r. n(t) is noise 
with a bandwith approximately equal to the receiver front-end 
bandwidth. Estimation of the delay (r) and Doppler (fo) is 
done by correlating the received signal with a model 

s m (t) = p(t — f)e _27r ^ c+ ^ jD ^ J (2) 


and determining f and fn which maximize this cross- 
correlation. 

This same correlation process can be applied to signals 
reflected from a random rough surface which are composed 
of an infinite distribution of reflected ray paths, each taking 
the form of (1) with a differential amplitude. The resulting 
expected value of this cross-correlation as a function of f and 
Id, 


Y (f,f D ) ) = (s(t)s* m (t-fJ D ) 


( 3 ) 



Fig. 1. Surface Scattering Coordinate Systems 


is the DDM, which can be shown to be a surface integral 
dependent on the PDF of surface slopes. The relevant geometry 
for this model is shown in figure 1 in which the GPS satellite 
is assumed to be sufficiently far away such that the signal 
incident on the surface can be described by a plane wave at 
the same elevation angle, 7, as the satellite. Figure 1 defines 
a coordinate system fixed in the scattering plane, the earth is 
assumed to be flat at the specular point, p is the position vector 
of a point on the surface, as measured in this plane. The slope 
at each reflecting point on the surface is thus given by the 
vector q(p) = ^-n — m, where rh and n are unit vectors in 
the incident and scattered directions, respectively. 

Using the geometric optics approximation, which requires 
that the radius of curvature of surface features are much larger 
than the GPS wavelength, the DDM model, derived first in [8] 
is given by: 


(|Y(r,/)| 2 ) 


2 f [ D 2 (p)a 0 (p) 
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Ti is the coherent integration time, D(p) is the footprint func- 
tion of the antenna, Rq, R are distances to the transmitter and 
receiver respectively, and x(Ar(p), A/(p)) is the ambiguity 
function of the transmitted GNSS signal (the correlation of (1) 
and (2)) where f s (p) is the doppler frequency at the scattering 
point. Dependence on the slope PDF enters through the bistatic 
scattering cross-section, cro(p). 
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where F is the complex Fresnel reflection coefficient and 
is the slope PDF. As shown in [9] this model may be rewritten 
as a convolution of A 2 with p(rj), the distribution of reflected 
power from the surface at delay 77. 


p(v) = Tf J J 


D(p)6(v -Ro-R) 
AnR 2 0 (p)R 2 (p) 


cr{p)d 2 p 


( 6 ) 


Surface slopes are generally assumed to follow a bivariate 
Gaussian distribution with extra terms to account for skewness 
and peakedness effects [13] [14] [15]. For this work we neglect 
the expanded terms and assume the slope distribution to be 


Gaussian and defined by the upwind and cross wind variances. 
Using a DDM for a fixed doppler allows retrieval of MSS 
without wind direction since the a one dimension slice of the 
DDM is not sensitive to the principal axis direction thus the 
total slope variance is estimated by assuming an isotropic slope 
distribution in the forward model. 

III. Experiment Description 

Experiments were conducted on 02 March 2009 flying 
on the NASA P-3 aircraft out of the Goose-Bay, Labrador 
airport to an area located at approximately 46 E, 57 N over 
the Atlantic Ocean. Data from three instruments, PALS [16], 
POLSCAT[17] and the DMR [9] [18] were used to record 
brightness temperature, wind speed from Ku-band polarimetric 
backscatter, and the GNSS-R DDM. A brief description is 
provided for each instrument and the flight campaign. Section 
IV will describe the data reduction and the inter-comparison 
of these measurements. 

A. Delay Mapping Receiver 

Since initial development in 1997, a number of functionally 
identical GPS Delay Mapping Receivers (DMR) have been 
built by the NASA-Langley Research Center and used in 
various field campaigns and studies. These include the NASA- 
U. S. Dept, of Agriculture Soil Moisture Experiments of 
2002 (SMEX02) [19][20][ and 2003 (SMEX03)[21][22] over 
agricultural areas, and flights over the desert environment 
Jornada Range near Las Cruces, New Mexico The DMR has 
also been used by the U. S. Army’s Cold Regions Research and 
Engineering Laboratory (CRREL) for surface snow and ice 
detection studies and has been flown onboard National Oceanic 
and Atmospheric Administration (NOAA) ’’Hurricane Hunter” 
aircraft for monitoring ocean surface wind speed [15] [23]. In 
the DMR, a bank of PRN code correlators is utilized to simul- 
taneously correlate the received surface-reflected signal over 
a range of 14 successive delay times (or delay bins). There is 
1/2 PRN code chip delay spacing between each correlator and 
the overall delay of the first correlator in the bank is set to 
ensure that the primary correlation peak (from signals near the 
specular point) is near the middle of the correlator bank (bin 
7). Following code correlation, further instrument processing 
produces a set of 14 ’power vs. [relative] delay’ measurements, 
recorded at 1/10 second intervals. Each correlation result is 
squared as part of instrument signal processing, and the set of 

squared correlations, Cj = ^ Y(fj,/^) ^ j = [1, ..14], is 
recorded. 

B. PALS 

The PALS instrument includes a combined L-band ra- 
diometer and scatterometer, operating at 1.413 GHz and 1.26 
GHz, respectively. It was built to investigate the benefits of 
combining passive and active microwave sensors for ocean 
salinity and soil moisture remote sensing. When it was built 
in 1999, it had one receiver at L-band and one receiver 
at S-band for either radar or radiometer. There were front- 
end microwave switches to sequentially select the antenna 


vertical and horizontal polarization ports to measure vertically 
and horizontally polarized radio emissions or radar echoes. 
To investigate the polarimetric observations of land or sea 
surfaces, the original S-band receiver was converted to another 
L-band receiver to enable simultaneous reception of vertical 
and horizontal polarizations. This conversion was completed in 
2006-2007 prior to PALS’ participation in the Cloud and Land 
Surface Interaction Campaign (CL ASIC) in 2007 [24] [16]. The 
PALS radiometer electronics were also modified to include a 
microwave hybrid after the front-end low noise amplifiers in 
the V and H receivers to produce the brightness temperature 
measurements at +45 and -45 degree linear polarizations, 
which can be used to derive the third Stokes parameter 
measurement. To achieve good calibration stability, the PALS 
front-end electronics were temperature-controlled to 0.1 C. We 
also applied temperature control (0.1 C) to the antenna with 
heater pads and temperature controller to ensure the stability 
of antenna gain. The PALS radiometer outputs include the ver- 
tically polarized (V) brightness temperature (TV), horizontally 
polarized (H) brightness temperature (T#) and the third Stokes 
parameter (U). The PALS radar provides the normalized radar 
backscatter cross-section (cr 0 ) for V-transmit/V-receive(VV), 
V-transmit/H-receive (HV), H-transmit/H-receive (HH), and 
H-transmit/V-receive (VH). 

C. POLSCAT 

The POLSCAT instrument is a Ku-band (13.9 GHz) 
scatterometer[25][17]. It includes two-axis gimbals for a 
conically scanning, parabolic antenna at constant elevation 
angle, which is controllable from 0 (nadir) to 65 degrees. 
POLSCAT transmits at a pulse-repetition-interval (PRI) of 
350 microseconds with alternating vertical (V) and horizon- 
tal (H) polarization from pulse to pulse. For each transmit 
pulse, POLSCAT detects the V and H polarization echoes 
simultaneously through two parallel receivers. This transmit 
and receiving design results in four polarization combinations, 
VV, HV, HH, and VH, for the detected radar echoes. The 
first letter of each polarization combination denotes the re- 
ceive polarization and the second letter denotes the transmit 
polarization. More detailed characteristics of the POLSCAT 
electronics, including the calibration loop for transmit power 
and receiver gain calibration, can be found in[25][17]. 

D. Data Collection Campaign 

Flights consisted of inbound and outbound paths and a 
variety of flight maneuvers including wing wags and constant 
roll turns. Aircraft altitude during the measurements was kept 
close to constant at approximately 3400 m. The flight path 
is shown in figure 2. On March 2 PolSC AT measured wind 
speeds ranged from approximately 4.5 to 29.5 m/s steadily 
increasing and decreasing over the inbound and outbound 
paths of the flight respectively. The DMR also recorded 
measurements of the DDM for subsequent retrieval of MSS, 
these results are shown in 4. Accurate aircraft orientation and 
positioning was provided by an Applanix GPS/INS system 
onboard the aircraft. This data allowed the full scattering 
geometry to be constructed for MSS estimation. The surface 




Fig. 2. Map of flight Path for March 2, 2009 


salinity was assumed to be constant over the flight areas and 
the expected impact of surface salinity and surface temperature 
was negligible over the range of brightness temperatures 
observed. 



Fig. 3. Correlation Data with best fit 


IV. Data Procesing 

A. Pre Processing 

The DDMs as measured by the DMR are the basic mea- 
surement used to retrieve surface slope statistics. The set 
of correlators used in the DMR, with fixed half-code chip 
spacing, is adjusted in instrument software to match the delay 
associated with the specular point. However, due to the fixed 
spacing changes in the path length of less than a half chip 
are not accounted for. Correcting the DDMs with the actual 
location of the specular point and using the known offset to 
the first correlator from the prompt correlator gives the relative 
distribution of power in the correct locations. An example of a 
re-aligned DDM is shown in figure 3, in this case the estimated 
MSS corresponding to the best fit was found to be .023. The 
details can be found in [5]. 

B. Surface Statistics Retrievals 

After aligning the correlation measurements, the power and 
corresponding delay is suitable for surface statistic retrieval. 
The estimation of MSS is done with a least squares approach to 
determine the slope standard deviation (cr), a delay offset (to), 
and a scale factor ( Sf act ). The model DDM, generated from 
the scattering model described above, is fit to the measured 
data, using a 3-element parameter vector p = {cr, To, Sf act }. 
Estimation of (to) and ( Sf act ) allow the model DDM to adjust 
for small uncertainties in delay as well as scaling to each 
individual block of measurements. 

P = argminj [Y 2 (r, f)-Y* (r, /; p)] T [Y 2 (r, /) -Y 2 U (r, /; p)} 

P 

(7) 

RMS slope (cr) was retrieved from each batch of 50 DDM’s, 
producing an independent estimate using 5 seconds of data. 
The DDM with the estimated waveform is shown in figure 3. 


V. Results 

Total MSS (cr 2 ) retrieved by performing the least squares 
estimation over the entire flight duration is shown in figure 
5 along with piecewise curve fitting solutions, the range of 



Fig. 4. cr 2 Estimate for March 2 2009 

cr 2 was approximately .0155 < a 2 < .0375. The solution of 
T\ 'fj to MSS was assumed to be of a form similar to that 
derived for wind speed [15], in this case a linear portion 
for small MSS and an exponential relation for large MSS 

Tb(Hy){? 2 ) = aeba + ceda with coefficients a, 6, c, d to be 

determined. 

In computing the exponential fits to the data certain points 
were excluded based on the cosine of the angle between the 
antenna boresight vector (b) and the scattering vector at the 
specular point (—h s ) as shown in figure 1. This was done 
due to possible nonuniform antenna gain pattern effects which 
are not taken into account in the forward model. This may 
also account for the varying groundtrack of the GNSS-R 
measurements relative to the aircraft fixed radiometer. This 
was computed by taking the dot product of the two vectors, 
the former of which was computed using the roll, pitch, and 
heading of the aircraft and is shown in figure 8. In figure 6 the 
residuals computed from an initial unweighted exponential fit 
to the data are plotted against the dot product value suggesting 
an empirical exclusion rule of (cos (6) < .851) to remove 
the largest errors. In figure 5 the least square curve fitting 
solution is shown where the influence of this rule can be 
seen, the points shown with a red x have been removed from 
consideration. The piecewise fits to MSS are given by 8 and 
9. 





Fig. 7. Tk vs. Tk stars represent a fused MSS and w po i estimation, dots are estimations using MSS alone LeftiT^y Right: 
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Fig. 5. Brightness Temperature with Piecewise Fit to cr 2 Top: Vertical 
Polarization Bottom horizontal Polarization 



Fig. 8. cos(9) on March 2, 2009 



40.4a 2 + 120.3 
213.8e" 249 1<T2 + lll.le 3 ' 6<T2 


for cr 2 < .0228 
for a 2 > .0228 


TbH 


107. Ict 2 + 72.9 

577.4e _319 - 2lT2 + 64.6e e - 5<j2 


( 8 ) 

for cr 2 < .0221 
for cr 2 > .0221 ( 


A multivariable solution was also found in both wind speed as 
measured by POLSCAT and MSS as measured by the DMR, 
a similar piecewise linear and exponential solution was used 
and the numerical coefficients in 10 and 1 1 were computed. 


Tbv 


. 99ct 2 + .18 w po i + 119.98 for cr 2 < .0228 
< .288 w poi + 33763. 4e 0012lT2 + 70.9e- 878 - 3<j2 
+3.78e- 11 ' 9<T2 ^° i - 33645.6.0 for cr 2 > .0228 

( 10 ) 


Fig. 6. Residuals vs. cos (6) Top '.T^y Bottom: 


Fit Type 

RMS Error [°K] 

Tbv ( w poi ) 

0.503131 

TbH ( w pol ) 

0.463798 

T bv (<r*) 

0.688453 

T bH (^) 

0.768167 

TbV ? w pol ) 

0.481308 

T bH ((T 2 ,Wp 0 l) 

0.462640 


TABLE I 

Summary of RMS Errors 


{ -17.9a 2 + .376 w pol + 72.9 for a 2 < .0221 
.388 w po i - 550097. 7e _866 1<T2 + 292.4 e - ° 51<j2 
+12.5e~' 04<72wpo ‘ - 232 for ct 2 > .0221 

( 11 ) 

In this work the antenna gain pattern was modeled as having 
a cos 2 dependence on the off nadir angle, D = (ph s ). In which 
b is along the antenna boresite and — fi s is in the directon 
of the specular point. Figure 6 shows the RMS error of the 
fit to Tty and T^h from MSS where there appears to be a 
dependency on the dot product between the antenna boresight 
vector (b) and the scattering vector at the specular point (—n s ). 
Figure 7 presents a scatter plot for PALS measured T b vs the 
estimated T ' b based on the curve fitting solutions to MSS (red) 
and the multivariable fit (blue) in both MSS and wind speed 
as measured by PolSCAT, for reference the 45° line is plotted 
in black. 

VI. Conclusions and Future Work 

For each type of fit for T b the resulting RMS error is shown 
in table I. Data from this experiment suggests that the shape of 
the delay-doppler map generated from GNSS-R signals may be 
used to estimate roughness effects on brightness temperature 
retrievals from microwave radiometry. Furthermore, current 
methods of correcting brightness temperature measurements 
due to surface roughness may benefit from coincident mea- 
surements of ocean surface statistics from reflected GNSS 
signals. There appears to be some dependency of RMS error 
on 0 indicating that a more accurate model for antenna gain 
or correction for different observation points between GNSS- 
R and radiometer groundtracks. Furthermore, using a full 2-D 
DDM to retrieve upwind and cross wind cr 2 and wind direction 
may improve the estimations. 
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